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We describe two peculiar galaxies falling into the massive galaxy clusters Abell 
1689 (z«0.18) and 2667 (z«0.23) respectively. Hubble Space Telescope images show 
-".^ , extraordinary trails composed of bright blue knots (-16. 5< M <-11.5 mag) and stellar 

• streams associated with each of these systems. Combining optical, near and mid- 

Oh' infrared and radio observations we prove that while both galaxies show similar ex- 

^ . tended trails of star-forming knots, their recent star formation histories are different. 

' One (« L*) is experiencing a strong burst of star formation, appearing as a rare ex- 

I ample of a luminous infrared cluster galaxy. In comparison, the other (w O.IL*) has 

Ci • recently ceased its star formation activity. Our model suggests that the morphologies 

and star formation in these galaxies have been influenced by the combined action of 
tidal interaction (likely with the cluster potential) and of ram pressure with the intr- 
' acluster medium. These results can be used to gain more insights to the origin of SOs, 

5-H ' dwarf and ultra-compact dwarf (UCD) cluster galaxies. 

! 
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1 INTRODUCTION 

Clusters are natural laboratories to study th e effects 
of en vironment on the evolution of galaxies (|Dresslerl 
Il980l ). A plethora of evidence shows that the proper- 
ties of late-type galaxies depend strongly on environ- 
ment: besides the well known morpholog y-density rela- 
tion (iDressleil llQSOl : IWhitmore et all Il993l ). in local clus- 
ters {z ^0.03) spiral galaxies are deficient in neutral 
hvdrogen ijCiovanelli fc HavnesI Il985l : ICavatte et al] Il990l : 
iBravo-Alfaro et al.ll200^ and have lower star formation ac- 
tivity than galaxies of similar t ype and size in low den- 
sity environments l|Lewis et al.l |2002| : I Gomez et al.l l2003l : 
iGavazzi et al.l |2006| ). Various physical mechanisms have 
been proposed to explain the different evolution of late 
type spirals in clusters and in the field. In general, they 
invoke either dynamical interactions of cl uster galaxies 
with the hot intracluster medium (ICM, iGunn fc GottI 



1 19721 : II arson et al.lll98(] | ), or grav itational interactions with 
nearby companions ("Mcrritt' 1983) , with the potential of the 
cluster ijBvrd fc Valtonen 199a:Vallurill993h . or with a com- 
bination of these two ( Moore et al.ll 19961 ). Interactions with 
the ICM are likely to be the dominant process at the present 
epoch and can account for the trunc ation of the gas disks in 
members of several local clusters (see lBoselli fc Gavazzill2006l 
and references therein). However, ram-pressu re cannot pro- 
duce the strong morph ology-density relation l|Dresslerlll980l : 
IWhitmor e et al.lll993l ). nor can it thicken the disk of a spi- 
ral galaxy and transform it into an SO (i.e. Hinz et aLll2003l : 
IChristlein fc Zabludofj|20oi : iBosehi et al.ll2006l ). 
This apparent contradiction could be solved if the struc- 
tures form hierarchically; that is, if galaxy clusters form not 
by accreting individual galaxies randomly from the field, 
but rather by the infall of less massive groups along fil- 
aments. These infalling groups have velocity dispersions 
that are much smaller than that of the cluster as a whole. 
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permitting the slow, st rong interactions normally associ- 
ated with field galaxies llFuiitall2004l : lMihosll2004l : iDresslerl 
I2OO4I : ICortese et al ] l2006bl ). Therefore a plausible evolu- 
tionary history would take into account that environmen- 
tal conditions and the physical properties of galaxies are 
changed significantly during cosmic time, changing the in- 
fluence o f various physical mech anisms on the evolution of 
galaxies (jBoselli fc Gavazzill2006l ). However, this hypothesis 
is far from confirmed, since we lack detailed understanding 
of the range of environmental e ffects that act as a func- 
tion of the age of the Universe (iDressleilliooi ). Ahhough 
star-forming galaxies in clusters at intermediate redshift ap - 
pear more disturbed (jOemler et al.|[l997l : Couch et al. 19981) 
and have h i gher star formation activity l|Butcher fc Oemlerl 
1 19781 . 11984 iFadda et al.ll200ol ) than local disk galaxies, it 
is still an open question which mechanisms are at play 
and how they influence the evolutionary history of cluster 
galax ies (|Balogh et al ] | 19991 : iDressTer' et al. l ll999l : lTreu et al.l 
|2003| ). To solve this riddle, we need to observe galaxies 
that physical circumstances and chance have revealed in 
rare moments of transformation. These peculiar systems 
can be used to probe different environmental effects and 
to constrain models of the evolutionary history of galaxies. 
Much of our kn owledge on the evolution of nearby galaxies 
in both groups (IDuc et aLll2OO0l : Ugl esias-Paramo fc Vflchej 



2001"; 'Sulentic et al."200l') and clusters 



Cortcsc et al. 2006b; Gavazzi et al.l 12001. : Kennev et al.l 
i995l : IVollmer et al.lf200ll . l2004h hasTn fact come from study- 
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ing of such systems. Unfortunately this information is dif- 
ficult to extend to high redshift because both clusters and 
;alaxies have evolved significantly: clusters were less relaxed 
Jehema et"aLll2005h and galaxies had higher gas content. 
Therefore, the effects of the same environmental mechanisms 
could depend strongly on the age of the Universe. 
In this paper, we present a multiwavelength analysis of two 
peculiar galaxies (hereafter referred as 235144-260358 and 
131124-012040) falling into the centers of the massive clus- 
ters AbeU 2667 (2 RiO.23) and Abell 1689 {z «0.18). Both 
these systems are associated with extended trails of bright 
blue knots and diffuse wisps and filaments of young stars, 
features ob served so far only in one other galaxy at simi- 
lar redshift l|Owen et al.ll2006l ). These two objects have been 
serendipitously detected by looking at the WFPC2 and ACS 
images of massive clusters at z~0.2. Th e sample of c l usters 
consist of the 10 clusters discussed in ISmith et al.l (|2005| ) 
plus A2667, A2390 and A1689. We therefore found 2 galax- 
ies with extended trails within 13 studied clusters all located 
at 0.175< z <0.25, suggesting that we are observing a very 
short snapshot of a critical stage in the evolution of these 
cluster galaxies. 

Because these two systems have significantly different opti- 
cal luminosities (« L* and «0.1L*) but are at similar dis- 
tances from the cores of clusters of similar mass, they rep- 
resent an interesting case for a comparison of the effects of 
similar environments on different-sized galaxies. 
Throughout this paper we assume a cosmology with Q,m = 
0.3, Q,\ = 0.7 and Ho = 70 km/s Mpc~^, implying a dis- 
tance modulus of 39.74 (40.33) mag and a linear scale of 
3.16 (3.71) kpc/arcsec for A1689 (A2667). 



2 THE DATA 

2.1 Optical Photometry 

The optical photometric data for this paper are extracted 
from deep HST images of AbeU 2667 and AbeU 1689. AbeU 
2667 was observed on 2001 October, using the Wide Field 
Planetary Camera (WFPC2) for total exposures of 1200 sec- 
onds through the F450W filter, and 4 000 seconds in F606W 
and F814W (see Fig. □ and FigEJ (|Covone et all l2006al ). 
The 3 a detection limit for point sources is « 26.00, 26.00 
and 25.00 mag in the F450W, F606W, and F814W bands, 
respectively. Deep observations of Abell 1689 were obtained 
from the ACS Guaranteed Time observations in 2002 June 
(see Fig. [land Fig|3]). A total of 20 orbits (Ril3.2 h) were 
taken in the three passb ands F475W, F625W and F850LP 
(|Broadhurst et ai]|2005l ). The 3 a detection limit for point 
sources is « 29.70, 29.20 and 28.30 mag in the fiUers F475W, 

F625W, F850LP res pectively. 

We used SExtractor (|Bertin fc Arnout j|l996l ) to detect and 
analyze the sources. For source detection, we used an image 
averaging the three band-passes; magnitudes were then de- 
termined from aperture photometry on the separate images 
for each filter. All magnitudes are in the VEGAMAG sys- 
tems. No correction for Galactic extinction was performed 
{Ay ^0.07 mag). Surface brightness profiles for these galax- 
ies were computed using the task ellipse within IRAF. The 
ellipticity and position angle were determi ned using the 
I ban d images following the procedure of ICavazzi et al.l 
(|2000l) . The disturbed morphologies of 235144-260358 in 
A2667 could possibly affect the shape of the surface bright- 
ness profiles at large radii (in particular at shorter wave- 
lengths), but not in the central regions where both the ob- 
jects still present a reasonably symmetrical shape. This is 
not the case of the edge-on spiral 131124-012040 which does 
not show strong asymmetries within the optical radius. 



2.2 Near Infrared Photometry 

Near-infrared H band observations for Abell 2667 and Abell 
1689 were obtained with ISAAC on the Very Large Telescope 
(VLT) in the spring of 2003 (ESO Programs 071.A-0428, 
PI. Kneib and 067.A-0095 PI. Knudsen), under photometric 
sky conditions with a mean seeing of «0.41 arcsec and ~0.58 
arcsec respectively. The total exposure time of 6529s for each 
cluster corresponds to a 5cr detection limit for point sources 
of ~24.6 ma g. All these observat ions have been reduced as 
described in lRichard et al.1 (|2006l ). 



2.3 Mid Infrared Photometry 

Spitzer imaging observations of Abell 2667 and Abell 1689 
were obtained as part of the GTO Lensing Cluster Survey 
(program 83, PI G. Rieke). MIPS (|Rieke et al.ll2004l ) 24/im 
images were obtained in photometry mode, with a total ex- 
posure time of «2700s. The data were proces sed and mo- 
saickcd following the procedures described in lEgami et al.l 
( 2006 ). Point source extrac tion and phot ometry were per- 
formed using DAOPH OT l|Stetsonl [l987l ) as described in 
iPapovich et ai] (|2004h . A PSF was constructed from the 
best-measured 30 point s ources in th e field; the Tiny Tim 
model of the 24/im PSF (|Kristlll995l ) was used to compute 
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the aperture correction. 131124-012040 in Abell 1689 is not 
detected in MIPS images. We derived a 3 a limit using a 
photometry aperture radius of 6 arcsec, a sky annulus be- 
tween 6 and 13 arcsec and an aperture correction of 1.698 
IRAC (iFazio et al.l |2004| ) four-bands (3.6, 4.5, 5.8 and 
8.0/xm) imaging was also obtained, with a total exposure 
time of 2400s per band for each cluster. Basic calibrated data 
were combined using a custom IDL mosaicking routine. For 
A2667 photometry was performed within apertures of ra- 
dius ~ 8.3 arcsec and no aperture corrections were applied 
because they are small with this large an extraction aper- 
ture. On the contrary in A1689 photometry was performed 
within a smaller aperture to avoid light contamination from 
nearby sources. We adopted a radius ~ 2.4 arcsec, sky an- 
nulus between 2.4 and 7.2 arcsec and aperture corrections 
of 1.213, 1.234, 1.379 and 1.584 at 3.6, 4.5, 5.8 and 8.0^m 
respectively. 

2.4 Radio continuum observations 

We obtained a 20 cm radio continuum measurement of 
Ril.4mJy for the galaxy in A2667 from the 1.4 GHz NVSS 
continuum survey. As this survey offers rather poor spatial 
resolution (~45 arcsec), we also constructed a 20cm map 
using higher resolution data from the NRAO VLA data- 
archive. Two different observations are available on A 2667 
at this frequency: 1) it was observed in October 2001 for 
3590 seconds in correlator mode 4, with a bandwidth of 25 
MHz, and using the CD-configuration (due to the low dech- 
nation); and 2.) a 3620 second observation was obtained in 
September 2002 with the same correlator mode and band- 
width, but in the BC configuration. We applied standard 
VLA calibration and imaging procedures with AlPS. We 
then combined the two data sets in the UV-plane using DB- 
CON. Images were generated with the task IMAGR and a 
weighting option robust=0, producing a map intermediate 
between natural and uniform weighting. After CLEANing, 
the resulting continuum map has a beam size of 16.7 x 13.1 
arcsec and an average rms noise of 0.12 mjy/beam. This 
map is shown in Fig[6] superposed on the HST image. 
The cluster A 1689 was observed with the VLA for a to- 
tal of 17405 seconds in the A-configuration in November 
2000 and March 2002. Combining the data sets using nat- 
ural weighting (to improve the sensitivity) we produced a 
map with a beam size of 2.1 x 1.6 arcsec, and an average 
rms of 0.15mJy/beam. No emission is detected at the posi- 
tion of the infalling galaxy. Taking a conservative detection 
threshold of 6a we estimate an upper limit of 0.90 mjy for 
the 20 cm radio continuum flux. 



2.5 Optical spectroscopy 

We obtained optical spectroscopy for 131124-012040 in Abell 
1689 as part of a wide fi eld (~ 30' x 3 0') spectroscopic sur- 
vey of the whole cluster l|Czoskell2004 Czoske et al. 2007, in 
preparation) using VIMOS on the VLT (Program 71.A-3065, 
P.l. Czoske). The LR-Blue grism was used, which provided a 
resolution of i? ~ 200 over a wavelength range from 3750 A 
to 6750 A. The dispersion was 5.35 A per pixel. We ob- 
tained three exposures of 840 seconds, for a total of 42 min- 
utes. The data were reduced using VIPGI jScodeggio et al.l 
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l2005h on site at lASF in Milano. The reduction involved 
bias subtraction, identification of the spectrum on the two- 
dimensional spectral image, interactive wavelength calibra- 
tion from observations of arc spectra and optimal extrac- 
tion using the method of iHorng l|l986l ). The spectrum has 
been flux-calibrated from observations of a spectrophoto- 
metric standard star, Feige 110. 

In addition we observed A 2667 and A1689 in June 2006 
with the LRIS instrument (|Oke et al.lll995l ') on Keck I. A 
5600 A dichroic separated the red channel of the instru- 
ment, equipped with a 400 1 mm"""^ grating blazed at 8500 
A, from the blue channel, equipped with a 600 1 mm~^ grism 
blazed at 4000 A. This setting covers the wavelength range 
3300 - 9200 A with a dispersion of 0.6/1.8 A and a res- 
olution of 4.3/6.3 A in the blue/red channel, respectively. 
On June 29, a 175"-long and l"-wide slit was aligned on 
the center of 235144-260358 in A2667 including two of its 
knots (Kl and K2, see Fig[SJ, using a position angle of 92.7 
East of North. Two exposures of 900 s were obtained un- 
der « 1.5"seeing. On June 30, a 30-slits mask was used on 
A1689 in order to target multiple-imaged candidates at the 
cluster center. One blue knot associated with the disrupted 
galaxy was included in a 9"-long and l"-wide slit from this 
mask (see Fig[S}. Four exposures of 1800 sec each have been 
obtained with an average seeing of 1.0". All these spectro- 
scopic data were reduced using standard IRAF procedures 
for flat-fielding, wavelength and flux calibrations, sky sub- 
traction and extraction of the spectra. 

Under visual inspection of the spectra we carried out the 
measurement of the emission and absorption lines using the 
task SPLOT in IRAF. For the spectrum of 235144-260358 
in Abell 2667 we de-blended the unde rlying absorpti o n from 
the H/3 emission lines as discussed in iGavazzi et al.l (|2004h . 
We evaluated the Balmer decrement from the ratio H/3/Hq 
(assuming r=10 000 K and n=100 e /cm^ . IOsterbrocklll989h 
and derived the correc ted line fluxes, relative to H/3, using 
the de-reddening law of lLegueux et al] ((1979.) . The observed 
Hq/H/3 ratio for 235144-260358 is «3.65 implying a gas at- 
tenuation A{Ha) «0.56 mag and a stellar continuum atten- 
uation j4(V|}_«0. 31 mag (assuming the Galactic extinction 
curve. [Peilll992l ). 



2.6 X-ray Imaging 

We downloaded the 9.2 ks Chandra ACIS observation of 
A2667 from the public archive and reduced it following the 
standard "threads" from CI AO data analysis software (Ver- 
sion 3.3j3- We searched the exposure-corrected images for 
sourc es, using the t a sk wa vdetect with angular wavlet scales 
from iBrandt et akl l|200ll ') and a signiflcance threshold of 
1 X 10~^. No source is detected at the position of 235144- 
260358. We measure an upper limit to the 2-8 keV flux of 
1.2 X 10"" erg s"^ cm'^ 



^ http:/ /cxc. harvard.edu/ciao/index. html 
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3 RESULTS 
3.1 Abell 2667 

Fig. E] (upper panel) shows an RGB image of 235144-260358 
in Abell 2667, and its global properties are summarized 
in Table [T] Its optical redshift is z ~ 0.2265, lying in 
the low velocity tail of the velocity distribution of Abell 
2667 (i. e. ~830 km s~'^ lo wer than the mean cluster 
velocity; ICovone et al.ll2006al ). This face-on galaxy lies at 
a projected distance of « 0.34 h^^ Mpc from the cluster 
center (assumed to coincide with the position of the central 
CD galaxy of A2667, see Fig. [TJ. This system is one o f 
the brightest galaxies in the cluster (jCovone et al.l [2006al ) ■ 
with both optical and near infrared (Mf450W «-21.50, 
Mh «-24.50) absolute magnitudes close to L* and a gas 
metallicitj0 of 12 + log{0/H) «9.0±0.1 (i.e. «1.4 solar 
metallicity). From the HST images, 235144-260358 appears 
to be a late-type galaxy (see Fig[3}, as confirmed by its 
structural parameters (see Table [l|. However this object 
is definitely not normal: it shows a disturbed morphology, 
with clear indications of stripping within its optical disk 
and a prominent one-armed spiral component as is typically 
obser ved in gravitationally perturbed systems (|Vollmerl 
l2003l ). Moreover, there is a significant nuclear enhancement 
in the optical surface brightness profiles (~ 2 mag within 
the central kpc: see Fig|4} , suggesting that it is experiencing 
a nuclear burst of star formation. This spike is particularly 
evident in the F450W band, where the central regions 
cannot be fitted with a simple deVaucouleurs profile. 
Spitzer observations of A2667 expand on the unusual 
properties of 235144-260358, since it is detected by both 
IRAC (3.6-8^m, see also FigO and MIPS (24 ^m). At 
z «0.23, the 8 ^m emission is dominated by a combination 
of the PAH bump (~ 6.2 /jm) and very small grains contin- 
uum l|Desert et al.lfl990l ). while the old stellar population 
dominates at shorter wavelengths. The observed 8/im/5.6 
Hm flux ratio ~ 6.3 (corresponding to the rest frame flux 
ratio 6.3/im/4.5 fj,m) i s consistent with the value observed in 
star forming galaxies (jPale et al.ll2005l ) , suggesting that the 
infrared emission is due to recent star formation activity. 
We used the X-ray data for a second test of whether the 
infrared emission is due to a burst of star formation or to an 
active nucleus (AGN). Comparing the X-ray upper limit to 
the 24 ^ m flux density (see Tab le [T]) with the help of figure 1 
of Alo nso-Herrero et al.l (|2004l ) confirms that this source is 
not AGN-dominated: its 2-10 keV/24/im fiux ratio is at 
least 4 times too low to lie within the range of typical AGN 
0. Finally a significant contribution from an AGN is ruled 
out by the emission line ratios obtained from the optical 
spectrum: log{[OIII]/H(3) ^-0.45, log{[NII]/Ha) « -0.33 
consisten t with the values typ ically observed in star forming 
galaxies (|Kewlev et ahlfioOlf ). We therefore used the 8 and 
24pm data to derive the total infrared luminosity, L{IR), 



^ The gas metallicity has been computed from the aver- 
age of five dif f erent empirical determ inations based on: iJ23 



Zaritskv etal] 11994 iMcGaughl Il99ll) . [ Nil] A6583/fOII ] A3727 
Kewlev fc Dopital I2OO2I) . [ NII]A6583/Ha jvan Zee et al] Il998t) 



and [OIII]A5007/[NII]A6583 llDutil fc Rovlll999h 
^ Neither o f these tests ca n exclude the presence of a Compton 
thick AGN HhTetaLllioO^ , but it is likely that the mid- infrared 
output of such objects is dominated by star formation. 



using the IR s pectra l ener gy distribution ( S EP) from the 
iDale fc Heloul l|2002l ) and icharv fc Elbaj ll200ll) librar y 
following the procedure described in lMarcillac et al.l l|2006l ). 
This method relies on the correlations between L(IR) vs. 
the luminosity at 7 ^ ,m and L{IR) vs. L(24/xm) shown in 
ICharv fc Elbaj (|200ll ). The SED templates were only used 
to interpolate at 8 and 24 fj,m. The resulting total infrared 
luminosity of 235144-260358 L{IR) Ri3(±0.25) lO" Lq, 
implies a current star format ion rate SFR ~ 53(±4.3) 
Mq yr~^ (using the relation of iKennicuttI 1 19981 ) . consistent 
with a SFR « 57 M© yr~^ ob tained from VL A continuum 
observations and the relation o f ICondonI l|l992l ). This galaxy 
is a rare example of a luminous infrared galaxy (LIRG) in 
a dense cluster. 

All the properties of 235144-260358 point out the pecu- 
liarity of its recent evolutionary history. However, the 
extended trails of bright blue knots, tracing its trajectory 
as it falls into the cluster core, make it truly extraordinary. 
A dozen such knots extend from the galaxy optical disk to 
a projected distance of ~ 80 h^Q kpc. Extended blue low 
surface brightness wisps and filaments lie along the same 
trail, supporting the hypothesis that all of these structures 
result from stripping (see Fig[3}. The knots have absolute 
F450W magnitudes in the range -16. 80< Mfasow <- 14.80 
mag, typical of d warf galaxies (jSandage et al.l Il985l ) and 
super star clusters (|Larsen fc Richtleilll999l ). and are barely 
resolved in the HST image implying an effective radius 
Te <0.45 h^Q kpc. 

The radio contours shown in Fig|6] appear elongated in the 
direction of the trail. A similar morphology seems also to 
be present in the Spitzer 8fim map shown in FigO which 
has the appearance of a head on the galaxy, with a tail 
tracing the current star formation associated with the blue 
knots. Moreover [OH] emission, not associated with any of 
the blue knots, extends from the galaxy for a total length 
of at least ~50 kpc (see Fig. [S]), suggesting the presence of 
diffuse ionized gas along the trails as already observed in 
nearby ram pressure stripped galaxies llGavazzi et ahl 1200 ll : 
lYoshida et al. I l2004l : ICortese et al.|[200m 
To constrain the ages of the blue knots, we compute the time 
evolution of the FA50W - F6QQW and ^6061^ - F814LP 
colors, using Starburst99 (|Leitherer et al.lll999l ). We assume 
a Salpeter IMF, solar metallicitjQ and two different star 
formation histories: an instantaneous burst and continuous 
star formation. For each star formation history we also 
compute a model including the contribution of strong emis- 
sion lines. We redshifted the synthetic Starburst99 spectra 
to the cluster distance and used the synthetic photometry 
package SYNPHOT in IRAF to compute the model colors 
in the WFPC2 passbands. In Fig. we plot a color-color 
diagram for the bright knots (black circles) to compare 
with the theoretical evolutionary tracks (solid and dashed 
lines). We present only the colors for the brightest knots: 
i.e. detected in all three HST bands (see Fig[3}. The arrow 
shows the effect of attenuation by dust on th e observed col- 
ors, a ssuming a Galactic attenuation curve (jSchlegel et al.l 
Il998l ). The models with emission lines (dashed lines) appear 



■* We also tested a Kroupa IMF and stellar metallicities in the 
range 0.004< Z <0.02, but the evolutionary paths do not signifi- 
cantly vary from the ones shown in Fig|9] 
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to fit the observed colors better than those without (solid 
lines). Most of the blue knots lie slightly below the modeled 
tracks for both an instantaneous burst and continuous star 
formation but are reasonably consistent with an age of the 
episode in the range 5 < t < 15 Myr in the first case and 
10 < t < 1000 Myr in the second. These values are probably 
upper limits since it is very likely that the star forming 
knots contain dust, as observed in e xtragalactic HII regions 
ijGerhard et al.l|2002l:[ Cortese et al ][2004) and star forming 



dwarf galaxies i Boselli et al.l 2003r Cortese et al.ll2006ah . 
Our spectra detect [Oil] in emission in both the knots (Kl 
and K2) included on the slit (see Fig. [8} confirming that 
these systems are still forming stars. For Kl we also de- 
tected [OIII] and Ha in emission {z ^ 0.227). The Ha flux is 
/ « 1.6 X 10"^^ erg cm"^ s"^ corresponding to a SFR ^0.02 
Mq yr~^ (not corrected for extinction). No continuum is 
detected above a flux limit / « 3 x 10~^^ erg cm~'^ s~^ A~^, 
implying a Ha equivalent width EW{Ha) ^ 50 A. The 
lower limit on EW{Ha) corresponds to an age of the knot s 
t ^5 Myr for an instantaneous burst (jLeitherer et al.lll999^ . 
This value is significantly shorter than the one obtained 
from the optical colors (see Fig. [Qjl implying a significant 
larger amount of dust in the knots {Av ~0.5-l mag) than 
the one observed in their parent galaxy {Ay «0.31 mag). 
In comparison, for a continuous star formation history the 
value of EW{Ha) corresponds to an age t ^1000 Myr, 
consistent with the estimate obtained from the optical 
colors, making this scenario much more likely. 



3.2 Abell 1689 

The disrupted galaxy in Abell 1689 is illustrated in Fig. 
[3] (lower panel) and its main properties are listed in Ta- 
ble 1. This galaxy lies at a projected distance « 0.24 h^^ 
Mpc from the cluster center (see Fig. [5)) and is ~ 2.5 mag 
fainter than the perturbe d galaxy in Abell 2 667 (i.e. with 
a luminosity of « 0.1 L*, Wilson et al] 1 19971 ). Its redshift 
is z RiO.1866 confirming that it belongs to A1689. Con- 
trary to 235144-260358, the surface brightness profile of 
this galaxy follows a typical exponential profile (see Fig|4}. 
However, the slopes of its color profiles are anomalous: in 
both F450W - FG25W and F625W - FSUW there is an 
inversion of the color gradients, with bluer colors toward 
the center. The galaxy outskirts have a F450W — FSUW 
color ~ 1.7 mag, ~0.6 mag redder than the galaxy cen- 
ter and consistent with the ty pical color of red se quence 
galaxies in the local Universe (|Bernardi et al.ll2003t ). Sim- 
ilar features have been observed in spiral galaxies in the 
Virgo cluster and su ggest recent (t g^30 Myr) gas strip- 
ping by ram pressure l|BoseUi et al.ll2006l '). 131124-012040 is 
neither detected at 24 fim by Spitzer nor in VLA contin- 
uum images (see Table [l]). This is consistent with the opti- 
cal spectrum of this galaxy (see FigO, which shows strong 
Balmer lines in absorption {EW{HS) ^6 A, D(4000) «1.21) 
and very little residual star formation {EW{[OII]) ~ 1.8 
A). This overall behavior suggests t hat the galaxy center 
has recently (t ^ 100 Myr, i.e . Po ggianti fc Barbarol 1 19971 : 
IShiova et al.ll2002l : iKauffmann et al]|2003l ) stopped forming 
stars. These spectral features are consiste nt with both a sim- 
ple truncated and a post-starburst SFH JShiova et al1l2002l : 
iPracv et al.ll2005l : ICrowl fc KennevlbOod T however the in- 



verted color gradients and the absence of a central enhance- 
ment in the surface brightness profile favor a ram pressure 
scenario (|Bekki et al.ll2005l : iBoseUi et al.ll2006l '). 
A « 30 h^Q kpc long trail, formed of at least six blue 
knots and a number of wisps and filaments, is associated 
with this system. The bright knots have absolute F475W 
magnitudes in the range -13. 5< Mf475w <-11.5 mag, ly- 
ing between dwarf galaxies and stellar clusters (« 3 mag 
fainter than the knots observed in Abell 266'i[3). The knots 
nearest to the galaxy are clearly resolved in the HST im- 
ages and have a typical size rs{F475W) RiO.8-0.9 kpc. In 
comparison, the most distant knots are not resolved im- 
plying a physical size re(-F'475W^) ^0.35 kpc. To determine 
the ages of the blue knots we computed the time evolution 
of the F475W - F625W and F625W - F850LP colors, as 
described in the previous section. The results of our anal- 
ysis are presented in Fig. |9] Most of the knots lie within 
the modeled tracks for an instantaneous burst with an age 
in the range 5< t <100 Myr, and are slightly above the 
model for continuous star formation with an age in the range 
10< t <1000 Myr. As for 235144-260358 no correlation is 
observed between the optical colors of the knots and their 
distance from the infalling galaxy (see Fig |10|l . The opti- 
cal spectrum obtained for the most distant knot (Knot A 
in Fig|8}, reveals the presence of strong [OH] in emission 
(/ ~ 5.4 X lO"'^^ erg cm~^ s"'^ A~^), while no continuum 
is detected at a limit « 5 x 10~^^ erg cm"'^ s~^ A~^, im- 
plying an EW[OII] ^ 108 A and showing that star forma- 
tion is still taking place in this system. Also Ha emission 
is detected, but it lies on a bright sky line and is affected 
by fringing making it impossible to use the Ha equivalent 
width to obtain an independent estimate of the age of the 
burst. It is interesting to note that the time scale necessary 
to invert the color gradients {t ^300 Myr) appears to be 
slightly longer than the age of the trails {t <100 Myr), sug- 
gesting that the two features could be signatures of different 
physical mechanisms. 



4 ENVIRONMENTAL EFFECTS ON THE 
EVOLUTION OF THE INFALLING 
GALAXIES 

The peculiar properties of the two galaxies falling into 
A2667 and A1689 suggest that both galaxies are undergoing 
strong transformations due to their interaction with the 
harsh cluster environment. However while these objects 
are at similar distances from the cluster centers and show 
similar extended trails of star-forming knots, their recent 
star formation histories are different. 235144-260358 is 
experiencing a strong burst of star formation, appearing 
as a rare example of a luminous infrared cluster galaxy. 
In comparison, 131124-012040 has recently {t < 100 Myr) 
ceased its star formation activity. To probe this difference, 
we investigate the effects of different environmental mecha- 
nisms on the properties and star formation history of these 
two galaxies. 

The high velocity dispersion of the two clusters 

5 ACS observation of A1689 are Ri3 mag deeper than WFPC2 
imaging of A2667. We cannot therefore exclude that knots as 
faint as the one in A1689 are present also in A2667. 
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(aiD >1000km ICovone et all l2006al : iLokas et all 

I2OO6I ) makes very unlikely a low velocity interaction or a 
merger with another cluster galaxy. This would not be the 
case if the two galaxies belong to smaller, kinematically 
distinct, dynamical units (e.g. infalling groups). However no 
observational evidences support this possibility. Therefore 
we will only consider high velocity galaxy-galaxy and 
galaxy-cluster gravitational interactions and ram pressure 
stripping by the hot intracluster medium (ICM) as possible 
mechanisms to explain the peculiarities of these two galax- 
ies. 

In order to reduce the number of free parameters in 
our model we assume that the two galaxies are falling on 
linear orbits into the cluster core. This very simple scenario, 
supported by the fact that infalling galaxies have usually 
highly eccentric radial orbits ( Bosclli & Gavazzi 2006), allow 
us to express the cluster-centric distance (r) as a function of 
the galaxy infalling velocity: 

(1) 



sin{arcos[y^i^)) 

where Vis is the (measured) velocity component along the 
line of sight and rproj is the cluster-centric distance pro- 
jected along the line of sight. Similarly, assuming that 
the trails of blue knots trace the galaxy's trajectory 
iMoore et al.lll999l '). their physical length is: 



(2) 



sin{arcos[y-^)) 
where Lproj is their projected length and their age is: 



^trails 



(3) 



Vinf all 

This value of ttraiis is based on the assumption that the 
trails are at rest with respect to the cluster, and must be 
considered as a lower limit for the real age of these features. 
Both cluste rs have a ID velocity dispersion aio 
1000km s"^ l|Covone et al.ll2006al : iLokas et al.ll2006l ). imply- 
ing a 3D infalling velocity Vinfaii ~ y/Saio ■ In the following 
we therefore assume a 3D infalling velocity in the range be- 
tween 1000< Vinfaii <1730 km s~^ (i.e. VSaio and am) 
that can be considered as upper and lower limit of the real 
value. The values so derived for the cluster-centric distance, 
the length and the age of the trails for the upper and lower 
limit of Vinfaii are summarized in Table [5] 

4.1 Gravitational interactions 

We can approximate the strengths of high velocity galaxy- 
galaxy and galaxy-cl uster interactions by u sing the im- 
pulse approximation (|Henriksen fc Bvrdlll996l l. The trans- 
verse and radial tidal accelerations experienced by the in- 
falling galaxy are: 



atr 



GMr,, 



R 



■[i?2 + (r + i?)2]i" 



Qrad — GMpert \_~2 



(4) 



(5) 



where Mp^rt is the mass of the perturber within r, i? is the 
radius of the perturbed galaxy (assumed to be «5 effective 



radii l|Gavazzi et al.ll2000l )) and r is its distance from the 
perturber. The radial tidal field tends to accelerate the edge 
of a galaxy. If it is more intense than the internal galaxy 
acceleration, given by 



a-gal — 



GMdyr, 



(6) 



where Mdyn is the dynamical galaxy mass, it is able to strip 
mater ial from the infalling galaxy. Following Gava zzi et al.l 
()1996h . we use the H-band rest frame luminosity of the two 
galaxies to estimate their dynamical masses within the opti- 
cal radius and derive their disk rotational velocities obtain- 
ing Mdyn ~ lO"'® M© and Mdyn « 10^° ® M© for 235144- 
260358 and 131124-012040 respectively. In the case of non- 
interpenetrating galaxy-galaxy interactions, the impact pa- 
rameter is at least equal or greater than the galactic radius 
(i.e. r ^ R) implying that material is stripped from an in- 
falling galaxy only if 



Mpert > 1.33 X Mdyn 



(7) 



(i.e. ^ 10"'^Mq and « 1O^°-^M0 for 235144-260358 and 
131124-012040 respectively). The possible perturber should 
not lie at a larger distance than the typical size of trails. In 
Abell 2667, the brightest objects within a projected distance 
of 100 kpc from 235144-260358 have an H band magnitude 
~ -24.5 mag (i.e. Mdyn ~ 10^^'^ Mq), fairly consistent with 
the lower limit required for effective stripping. Unfortunately 
their recessional velocities are unknown, making impossible 
a more detailed analysis of their possible interaction with 
235144-260358. In AbeU1689, the giant face-on barred spi- 
ral projected at ^20 kpc NE from 131124-012040 (see Fig. 
[2| is the only object (within 100 kpc) satisfying Equation [T] 
However the trail of blue knots is pointing in the opposite 
direction than the one expected in the case of interaction 
between the two objects (i.e. towards the perturber) and 
the galaxy has a redshift of z «0.1924 i.e. 1680 km/s higher 
than the recessional velocity of 131124-012040, making un- 
likely an interaction between the two objects. 
To quantify the effect of tidal forces from the cluster poten- 
tial well on an infal ling galaxy, we assume a NFW profile 
(|Navarro et al.|[l997l ) for the cluster mass distribution: 



M(< r) = Mo[lnil + -) - -^^^ 



where 



Mo = 47r|^ [nM{l + zf + Qa] 



for r < raC 



3^3 



200c'' r 



8ttG 



3[/n(l + c)-c/(l + c)] 



(8) 



(9) 



and r-g and c are the scale radius and concentration parame- 
ter of the mass distribution. The values adopted for the two 
clusters are summarized in Table [3] As shown in Fig llll for 
both our galaxies the radial acceleration from the cluster po- 
tential is higher than the internal acceleration for a cluster- 
centric distance smaller than « 0.45 h^g Mpc. Therefore, 
depending on their real infalling velocity, the two objects 
are at the edge or are just entered the region where material 
can be efficiently stripped by gravitational interactions. This 
simple calculation gives a lower limit for the real efficiency 
of mass loss, since higher rates will occur in the presence 
of substruct ures and infalling groups, as is likely in these 
two clusters (ICovone et al.ll2006bl: iLokas et al.]|2006l ). More- 
over, tidal heating (|Tavlor fc Babulll2001^ produced by the 
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varying clust er gravitationa l field will significantly acceler- 
ate mass loss (|Gnedinll2003al lbh. although it is not considered 
in our model. 

In contrast to the radial acceleration, which tends to strip 
material from an infalling galaxy, the transverse field com- 
presses the interstellar medium to produce gas infall to- 
ward the center and may trigger a burst of star forma- 
tion. Gas clouds experience a velocity perturbation due to 
the transversal tidal acceleration and collide with other gas 
clouds. The increase in the cloud velocity can be estimated 



V 



(10) 



where A(t) is the age of the interaction. The velocity in- 
crease and cloud collision produce a density enhancement in 
the centre of the galaxy, which is proportional to the Mach 
number (MN) squared. Consequently the critical mass for 
the cloud collapse (which is proportional to pgSs^) decreases 
by a factor MN~^ and, in the case of a strong perturbation, 
could become smaller than the typic a l mass of a galactic dis k 
HI cloud («300 M©, ISpitzej[l97i : |jog fc Solomoiil Il992l) 
favor ing new episodes of star formation ( Henriksen fc BvrdI 
1 19961 ). Fig.[T2]shows the ratio between the typical mass of HI 
clouds and the critical mass for cloud collapse in the case 
of high velocity interactions, assuming 

1O"-^M0 for 235144-260358 and 131124- 



and M 



pert 

012040 respectively as discussed above. It appears that in 
both cases high velocity galaxy-galaxy interactions are not 
enough strong to trigger a burst of star formation as ob- 
served in 235144-260358. 

This is not the case of galaxy-cluster interactions. Fig. 1131 
shows again the ratio between the typical mass of HI clouds 
and the critical mass for cloud collapse as a function of the 
cluster-centric distance in case of interaction with the clus- 
ter potential. In this case two galaxies are in two different 
regimes, whatever are the initial conditions in our model. 
While in 235144-260358 the critical mass is already below 
~300 M© and the compressed gas is able to collapse and pro- 
duce new stars, in 131124-012040 this is still not the case. 
This result is consistent with our observations and indicates 
that tidal forces from A2667 may have triggered the strong 
starburst in 235144-260358. 

In summary our model suggests that gravitational interac- 
tions with the cluster potential alone are able to strip ma- 
terial from the two infalling galaxies and to trigger a burst 
of star formation in 235144-260358. Even if we cannot com- 
pletely exclude a role of high velocity galaxy-galaxy inter- 
actions on the evolution of these systems, it appears clear 
that they cannot account for all the properties of the two 
infalling galaxies. 



4.2 Ram pressure stripping 

Although the tidal interaction hypothesis is consistent with 
the presence of a strong starburst only in 235144-260358, it 
is not able to explain why 131124-012040 shows clear signs 
of a recent truncation of its star formation. Both A1689 and 
A2667 are X-ray bright clusters suggesting that the effects of 
the hot intracluster medium could be significant. Therefore, 
to estimate the effects of ram pressure stripping on the in- 



falling galaxies, we adopt the classical iGunn fc GottI (|l972l ) 
criterion: 



Pram = pICMV ^ 2TTG'Sstar'Sgas 



(11) 



where picM is the density of the cluster medium, E^tar and 
Ejas are the galaxy stellar and gas density, and v is the 3D 
infalling velocity of the galaxy (here assumed to be in the 
range 1000< v < 1730 km s~^ as discussed in the previous 
sections). We use a /3 model density profile for the ICM: 



P^r) = pa 



[1 + (r/r2)]3/3/2 



(12) 



(the values adopted for the different clusters are listed in 
Table |3]). We assume that the stellar and gas distributions 
of our galaxies are exponential, as confirmed by their struc- 
tural parameters . The gas and stellar density profiles are 
(Domainko ot aLl l2006h 



^star 



star , gas 
'^'^-^Ostar,gas 



■ exp{ — r / Rostar.gas) 



(13) 



where Ro is the scale length of the exponential profile (i.e. 
0.59 Te). Assuming a gas scale length Rogas ~ 1-8 Rostar 
iGavatte et al.|[l994 ). a Mgas/ Mstar ratio » 1 , typi cally ob- 
served in late type galaxies (iBoselli et al.' '20 0l|), and a 
Mstar/ Lh ratio « 1 l|Gavazzi et al. 1996; Mc Gaugh et al.l 
l2000h the typical stripping radius is given by the following 
relation: 



Rstr.p ~ 0.64_Ro X In i 



G{Lh/Lc, 



1.8^ pic MV^2nR^^t^,.' 



(14) 



and the mass of gas stripped by ram pressure is: 



-I- l)exp{~Rstr,p/Ro) 



(15) 



M — f Rstrip 

Lq Ho 

In Fig. [14] (left) we show the variation of the stripping radius 
as a function of the distance from the cluster center, for three 
different values of the 3D infalling velocities assumed in our 
model {^lOOOkm s~\ «1410fcm s"^ and ^1730km s"^). 
While at its current location 131124-012040 has almost been 
totally stripped by ram pressure {Rstrip ^ 0.9re), in 235144- 
260358 ram pressure has only affected the outer galaxy re- 
gions {Rstrip ^ 1.2rc). The same result can be analyzed in 
terms of HI deficiencjQ (see Fig. [Til right). 131124-012040 
has already lost ^80 % of its original gas content and if ob- 
served in a local cluster it would be classified as a highly HI- 
deficient object. Conversely, ram pressure has only stripped 
a tiny fraction of the gas from 235144-260358, whose HI 
deficiency (~0.25) would approximately lie at the edge be- 
tween normal and deficient galaxies (HI-deficiency~ 0.2). 
We remark that the Hl-deficiency shown in Fig |14l is not 
determined from observations, but is obtained from our an- 
alytical model. 

Comparing Fig |14l to Fig llll it appears clear that ram pres- 
sure in 131124-012040 has become efficient before gravi- 
tational interactions were able to strip material from the 
galaxy. This is qualitatively consistent with the different 



^ The HI deficiency is defined as the diflerence, in logarithmic 
units, between the observed HI mass and the value expected from 
an isolated galaxy with the same morphological type T and op- 
tical linear diameter D: HI DEF = < log A/j^ / (T°'^ ) > 
-logM'^l teavnes fc Giovanelli[[l983) 
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time-scales of the interaction determined from the inversion 
of the color gradients {t ^ 300 Myr, likely produced by ram 
pressure stripping) and from the blue star-forming trails {t < 
100 Myr, clear signature of gravitational interactions). 

4.3 The origin of the blue star forming knots 

The mutual effects of gravitational interactions and ram 
pressure have already been observed in several cluster galax- 
ies, however the tails of blue star forming knots here dis- 
covered represent an extremely rare feature, to our knowl- 
edge previously observed only one other starburst gal axy in 
the AbeU cluster 2125 {z RiO.247, lOwen et al.]|2006l ). The 
morphology and luminosity of the knots suggest that we 
are dealing with dwarf gala xies and/or stellar super-clusters 
l|Fellhauer fc Kroupa|[200^ ) . This could explain the observed 
difference between luminosity of the knots («2.5 mag) in the 
two galaxies since the luminosity of the brightest star clus- 
ters is usually correla ted with the SFR of the parent galaxy 
l|Weidner et al.ll2004l ). 

The properties of the knots (i.e. colors and emission lines) 
suggest that they are undergoing an extended period of star 
formation as discussed in Section 3. From the model de- 
scribed in the previous sections the dynamical age of the 
trails results 50< ttraii < 150 Myr and 20< Uraii < 60 Myr 
for 235144-260358 and 131124-012040 respectively, fairly 
consistenlQ with the age inferred from their optical colors 
(see Fig.l^J. It is impossible to determine whether the knots 
where already forming stars when they have been stripped or 
if their activity was triggered by an external mechanism once 
in the ICM. We can exclude that these systems have been 
formed in the ICM by the accretion of unbounded material 
stripped by the parent galaxies: the combined effects of the 
cluster tidal field and ram pressure tend to inhibit the forma- 
tion of bound systems from the collapse of stripped material 
l)Mihosll200.j ). This scenario is consistent with the numeri- 
cal simulations. lElmeereen et all (|l993h showed that gravi- 
tational interactions can lead to the formation and ejection 
of peripheral self gravitationally bound clouds with masses 
^ 10* M0, which be gin their life in a majo r burst of star 
formation. Moreover iBekki fc CouchI (|2003h have recently 
demonstrated that ram pressure can trigger the collapse of 
stripped clouds leading to a burst of star formation, sug- 
gesting that even the formation and evolution of the blue 
star forming knots is probably driven by the mutual effects 
of gravitational interactions and ram pressure. Only deeper 
spectroscopic observations will shed light on the star forma- 
tion history of these rare objects. 



5 DISCUSSION & CONCLUSION 

The analysis in this paper allows us to propose a scenario for 
the evolution of the two disturbed galaxies in AbeU 2667 and 
AbeU 1689. These objects are currently falling into massive, 
gas-rich galaxy clusters with similar mass and gas density 
profiles (see Table |3]). Under the combined action of tidal 

However the value of ti^ail obtained in the previous section 
must be considered as a lower limit since it assumes that the 
stripped material is at rest with regard to the cluster. 



forces (more likely from the cluster potential) and of ram 
pressure by the ICM, their morphologies and star forma- 
tion are strongly perturbed. Self gravitational bound sys- 
tems are ejected from the main galaxies and stars and ion- 
ized gas are stripped from the stellar disks producing the 
observed tails of blue knots and stellar wisps tracing the 
infalling trajectory of these systems into the cluster core. 
Only the tidal field of AbeU 2667 is able to drive a gas in- 
fall into the center of 235144-260358 triggering a nuclear 
burst of star formation and making this galaxy a rare ex- 
ample of luminous infrared cluster galaxy. Simultaneously 
ram pressure stripping by the hot intracluster medium strips 
the neutral hydrogen from the galaxy outskirts but it is 
not able to affect the central regions where the starburst 
is taking place. Conversely in 131124-012040 gravitational 
forces are not strong enough to trigger the collapse of gas 
clouds while ram pressure is already extremely efficient. At 
the present galaxy location ram pressure has stripped at 
least the «80% of the original neutral hydrogen content, 
quenching the star formation activity in this object, as con- 
firmed by the strong B aim er lin e in absorption observed in 
the optical spectrum (jShiova et al.i 120021 ) and by the inver- 
sion of the optical color gradients along all t he galaxy extent 
(jBoseUi fc GavazzillioollBoseUi et al.il2006l l. 
A larger statistical sample is necessary to determine whether 
we are witnessing a common snapshot in the evolution of 
cluster galaxies or an extremely rare phenomenon. In fact, as 
discussed in the Introduction, only these 2 galaxies out of 13 
different clusters imaged at 0.175< z <0.25 show extended 
trails of blue knots. Within the WFPC2 field of view (RiO.25 
Mpc^ at 2 « 0.2) there are typically ~50 clu ster members 
but o nly the «20% of them are spiral galaxies (jBalogh et al.l 
l2002h . implying a frequency of «1.5% (2 over 130) among 
cluster spirals at z «0.2. This value is fairly consistent with 
the expected frequency roughly obtained dividing the typ- 
ical time scale of the interaction (^200 Myr) to the age of 
the cluster («11 Gyr). 

If we are witnessing a common step of cluster galaxy evolu- 
tion what can we learn by studying these two rare objects? 
AbeU 2667 and AbeU 1689 have comparable mass and gas 
profiles and the two galaxies are approximately at the same 
projected cluster-centric distance, suggesting that the abso- 
lute intensity of the cluster tidal field and of ram pressure 
by the cluster ICM is approximately the same in the two 
environments. Therefore we can speculate that the different 
recent evolutionary history of the infalling systems could be 
in part related to their different properties (i.e. their differ- 
ent luminosities: ~ L* and «0.1L*). In this case, our result 
suggests that giant spiral galaxies infalling into the core of 
massive clusters are mainly perturbed by the gravitational 
interaction with the cluster. Stars respond by forming arms 
and bars, while the gas flows directly toward the central re- 
gion within t «100 Myr. The sinking of the gas towards the 
center triggers a burst of st ar formation an d is able to al- 
ter the galaxy morphology (lono et al.ll2004l l. Ram pressure 
stripping produces a truncation of the disk but only in the 
outskirts of the galaxy being not efficient within the opti- 
cal effective radius. When all the remaining fuel has been 
consumed by the star formation this galaxy will not longer 
appear as a disky gas rich systems but more likely as a bulge 
dominated quiescent spiral. This is not the case for less mas- 
sive galaxies. Ram pressure is much more efficient on low 
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mass systems and it is able to strip a considerable fraction 
of the neutral hydrogen from the inner part of these galaxies, 
preventing the gas from sinking toward the center driven by 
tidal interaction and quenching the star formation history. 
Within Ril Gyr this object will not appear as a blue spiral 
any more bu t will probably look like an early type (e.g. red) 
disky spiral (jShiova et al.ll2()02 ). The different evolutionary 
scenarios for the evolution of low and high mass infalling 
galaxies emerging from our analysis apparently fit with re- 
cent observations and models, suggesting that the bulk of 
the cluster population of giant bulge dominated early type 
spiral galaxies can o nly be formed during some kind of gravi- 
tational interaction ('Dressier 20041: lMihosl l2004^ . while lower 
mass systems ca n be transformed by si mple gas removal from 
healthy spirals (|Poggianti et al.ll2004l ). 

The properties of the blue knots stripped from the in- 
falling galaxies deserve particular attention. These systems 
have a luminosity (-16.5^; M ^-11.5) and a physical size 
(re ^ 0.45 kpc) typical of dwarf galaxies and consistent with 
the ultra compact dw arf galaxies (UCD. lHilker et al.lll999l : 
[Phillipps ot al.' 2001), recently discovered in Abell 1689 
(Mieskc ct al. 2004, 200^). There are two comp eting forma- 
tion scenarios to explain the origin of UCDs. iBekki et al.l 
propose that they are the remnants of stripped dwarf 
galaxies. In this scenario a nucleated dwarf looses its en- 
velope and a great part of its dark matter content due 
to tidal interaction with a nother object. On the contrary 
iFellhauer fc Kroupal l|2002l ) propose that UCDa could origin 
from the amalgamation of rich aggregates of young massive 
star clusters that can form during gravitational interactions 
between gas-rich galaxies. It appears clear that the knots 
discovered here strongly support the second scenario, sug- 
gesting that at least part of the population of ultra-compact 
dwarfs originate from young massive star clusters: we are 
probably for the first time witnessing the dawn of the UCDs. 
This scenario is also consistent with the recent discovery of 
a massive extra galactic star cluster (M > 10® M©, t ~ 700 
Myr) lying at a projected distance of 17 kpc from the merger 
remna nt NGC3310 and likely formed during the merging 
event (iKnapp et al.ll2006l ). 

Finally the diffuse stellar streams and ionized gas observed 
along the trails suggest that the mechanisms acting here will 
significantly infiuence the properties of the intracluster light 
and contribute to the enrichment of the ICM. The results 
here obtained might be representative only of the clusters at 
z ^0.2 where the infalling rate is higher and galaxies have 
an higher gas content than the one observed in local clusters 
of galaxies. 
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Table 1. Properties of the two disrupted galaxies 





A1689 




A2667 




131124 - 012040 




235144 - 260358 


a (J. 2000) 


13:11:24.86 


a (J. 2000) 


23:51:44.03 


S (J. 2000) 


-01:20:39.9 


5 (J. 2000 


-26:03:59.6 


z 


0.1870 


z 


0.2265 


F475W (mag) 


-19.09 


F450W (mag) 


-21.54 


F625W (mag) 


-19.95 


FmeW (mag) 


-22.20 


F850LP (mag) 


-20.49 


F8UW (mag) 


-22.87 


H (mag) 


-21.82 


H (mag) 


-24.47 


F(3.6fim) (mjy) 


0.026 


F(3.6fim) (mJy) 


0.41 


F(4.5/im) (mJy) 


0.019 


F(4.5/im) (mJy) 


0.43 


F(5.8/im) (mJy) 


0.005 


F{5.8fim) (mJy) 


0.29 


F{8fim) (mJy) 


0.008 


F{8^lm) (mJy) 


1.81 


F(24/im) (mJy) 


<0.03 


F(24/im) (mJy) 


4.20 


F(20cni) (mJy) 


<0.90 


F(20cm) (mJy) 


1.40 


re (F850LP) {h-^ kpc) 


3.1 


re (F814W) {h-^ kpc) 


5.5 


fie (F850LP) (mag arcsec-^) 


20.07 


He (F814W) (mag arcsec"^) 


19.99 


C31 (F850LP) 


3.13 


C31 (F814W) 


3.02 



Table 2. Cluster centric distances and properties of the trails as a function of the 3D infalling velocity assumed in the model described 
in Section 4. 



Cluster 




fall «1730 km s ^ 




fall «1000 km s-i 




D 


^trails 


^trails 


D 


^trails 


^trails 




kpc 


kpc 


Myr 


kpc 


kpc 


Myr 


A1689 


280 


35 


20 


450 


60 


60 


A2667 


390 


90 


54 


610 


140 


150 



Table 3. Best fitting parameters for the total mass and gas density profiles of A1689 and A2667 



Cluster 


Mass Profile 




Gas Profile 






rs 


c Ref. 


P 


Rc 


po 


Ref. 




Mpc 






Mpc 


g cm^"' 




A1689 


0.31 


8.2 1 


0.55 


0.07 


1.77x10-24 


3 


A2667 


0.49 


6 2 


0.52 


0.049 


2.17x10-24 


3 



References: l lBroadhurst et al.l l l2005h 2) ICovone et all (Hooii); 3) lOta fc Mitsudal l l2004h . 
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Figure 1. HST-WFPC2 F450W image of Abell 2667. The position of 235144-260358 is indicated. 
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Figure 3. RGB color image of 235144-260358 in Abell 2667 (upper panel) and of 131124-012040 in Abell 1689 (lower panel). The white 
circles indicate the knots studied in Sec. 3 
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Figure 4. The surface brightness and color profiles in the three HST bands for 235144-260358 in Abell 2667 (left) 131124^012040 in 
Abell 1689 (right). 
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Figure 6. The VLA 20cm radio contours (2,4,8, 16, 34 and 64 a) of Abell 2667 superposed on tlie HST F450W image. 
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Figure 7. The optical spectrum of 235144-260358 in Abell 2667 (upper) and 131124-012040 in Abell 1689 (lower). The main emission 
and absorption lines are indicated. 




Figure 8. Left: The 2D spectrum of 234144-260358 in Abell 2667 and of two blue knots near the wavelength of the [Oil] emission line 

(lower panel). The slit position is indicated by the two solid lines in the upper panel. [Oil] diffuse emission appears clearly between the 
galaxy and the blue knots. Right: The 2D spectrum of knot A in Abell 1689 near the wavelength of the [Oil] emission line (lower panel). 
The position of the knot is indicated by the black circle in the upper panel. 
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Figure 9. Colour-Colour diagram for the blue knots in Abell 2667 (left) and Abell 1689 (right). The black circles with their associated 
error bars represent the knot colors and their uncertainties. The green squares indicate the colours of the two disrupted galaxies. The 
dashed lines are the model predictions for the evolution of SSPs (red) and continuous (blue) star formation accounting for the contribution 
of emission lines. The solid lines represent the same models without the contribution of emission lines. The numbers indicate the age of 
each model in Myr. The black arrows show the effect of internal extinction on the observed colours. 
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Figure 10. Left; Distribution of the F450W-F606W colours for the star forming knots in Abell 2667 as a function of the distance from 
235144-260358. The red points represent the average values for the nearest and most distant star forming knots. Right; The same diagram 
for the F475W-F650W color of the star forming knots associated with 131124-012040 in Abell 1689. 
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Figure 11. The ratio of the radial acceleration due to the cluster potential to the internal galaxy acceleration as a function of the cluster- 
centric distance for A1689 (upper panel) and A2667 (lower panel). The shaded regions indicate the range of cluster-centric distances 
assumed in the model (see Section 4). 
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Figure 12. The ratio of the typical HI cloud mass in galactic disks to the critical mass for gas collapse in the case of galaxy galaxy 
interactions as a function of the distance from the pcrturbcr for 235144-260358 in AbcU 2667 (solid line) and of 131124-012040 in Abell 
1689 (dotted line). The black circles indicate the values for the closest interaction possible (rsi Rgai)- The horizontal dotted-dashed line 
shows the transition between the stable {M^iay^^^/Mcritic <1) and unstable {M^i^y^^g/Mcritic >1) regime. 
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Figure 13. Same as Fig |12l in case of a cluster-galaxy interaction in A1689 (upper panel) and A2667 (lower panel), for three different 
values of the 3D infalling velocity: 1000 km s~^ (solid line), 1410 km s~^ (dotted line) and 1730 km s~^ (dashed line). The black circles 
indicate the current galaxies positions assumed in each model. 
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Figure 14. Left: The variation of the stripped radius as a function of the cluster-centric distance for 235144-260358 in Abell 2667 (lower 
panel) and of 131124-012040 in Abell 1689 (upper panel) for three different values of the 3D infalling velocity. Symbols as in Fig |13l 
Right: The variation of the HI deficiency as a function of the cluster-centric distance for 235144-260358 in Abell 2667 (lower panel) and 
of 131124-012040 in Abell 1689 (upper panel) for three different values of the 3D infalling velocity. The shaded area indicates the typical 
range of Hl-deficiency observed in unperturbed galaxies. Symbols as in Fig |13l For Vi„faii 1200 km the galaxy has lost all its 
original HI content. 



